Muon's Behaviors under Bremsstrahlung with both 
the LPM effect and the Ter-Mikaelian effect and 
Direct Pair Production with the LPM effect 

S. Polityko\ N. Takahashi^, M. Kato^ Y. Yamada^ and A. Misaki^ 
^Irkutsk State University, Irkutsk, 664003 Russia 
"^Facutly of Science and Technology, Hirosaki University, 

Hirosaki 036-8561 Japan 
^Kyowa Interface Science Co. Ltd., Asaka, 351-0033 Japan 
'^Hitachi Goverment and Public Cooperation System 
Engineering Ltd., Koto-ku, Tokyo, 136-8633 Japan 
^National Craduate Institute for Policy Studies, 
2-2 Wakamatsu-cho, Shinjyuko-ku, 162-8677 
and Institute for Policy Science, Saitama University, 
Urawa 338-8570 Japan 

November 1999 

Abstract 

Differential and integral cross sections of the muon are calculated in the materials: water, standard 
rock, iron and lead with and without the LPM effect. The corresponding cross sections are also 
calculated with dielectric supression effect (Ter-Mikaelian effect), in addition to the LPM effect. In 
our calculations the LPM effect for muon is provided to be effective from 10^** eV to 10^° eV depending 
on materials in the bremsstrahlung process, while it is provided to be completely negligible in the 
direct pair production process up to 10^^ eV even for lead. 

As for the dielectric suppression effect it is only effective between 10^^ eV and 10^'^ eV in the 
case with the LPM effect, while it is effective above 10^'^ eV in the case without LPM effect. To 
demonstrate the importance of the LPM effect in the bremsstrahlung process, depth intensity relation 
of muon, energy spectrum, range distribution and survival probability are calculated taking into 
account bremsstrahlung spectrum. 



1 Introduction 



In the early 1950's , Landau , Pomeranchuk and Migdal realized that because of the low longitudinal 
momentum transfer between the nucleus and the fast particle, bremsstrahlung is not instantaneous, but 
occurs over a finite formation zone. During this time, external influences can perturb the fast particle 
and suppress the photon emission. When this happens, the usual Bethe-Heitler formula fails. Initially, 
Landau and Pomeranchuk studied suppression by multiple scattering with semiclassical argument 
Later, Migdal presented a fully quantum treatment [Q. Now this phenomenon is well known as the 
LPM effect. Cosmic ray physicists have been interested in the LPM effect which becomes effective at 
extremely high energies. Besides the LPM effect is expected to change drastically the structure of the 
electromagnetic cascade showers at extremely high energies, the structures of which are closely related 
to their physical interpretation. Up to now, the experimental proofs of the LPM effect have been tried 
by many investigators in both cosmic rays |^ and accelerators 

Recently, the existence of the LPM effect has been finally confirmed by the SLAC electron accelerator 
and the longitudinal effect Q has been also reconfirmed . Electromagnetic cascade shower with the 
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LPM effect is called the LPM shower, while the LPM effect, namely that at not-extremely high energies 
is called the BH (Bethe-Heitler) shower. The LPM shower has two distinguished characteristics over the 
BH shower: a) The average behaviors of the LPM shower are quite different from those of the BH shower 
l^and b) The individual behaviors of the LPM shower are quite average behaviors of the LPM showers 

From the first characteristics of the LPM shower, we could easily understand that the size of the 
experimental apparatus for extremely high energy phenomena required should be much larger than that 
in usually imagined, while from the second characteristics of the LPM shower we understand that we 
should create a new method to analyze extremely high energy phenomena which could not be analyzed 
by the average theory. 

Up to now, we have discussed the characteristics of the LPM effect on the electrons and photons. On 
the other hand, the LPM effect due to muon is also important for extremely high energy phenomena., 
particularly, for studies on extremely high energy neutrino astrophysics which are quite an uncultivated 
field. Now, there are two test pilot experiments for high energy neutrino astrophysics experiment in Baikal 
12 . In future, we expect to construct the gigantic extremely high energy neutrino 
astrophysics detector, where the principle of the construction for the apparatus should be different from 
that of extremely high energy neutrino astrophysics due to the presence of the LPM effect. For the 
moment, it is more plausible that extremely high energy neutrino could be detected as extremely high 
energy muon from the neutrino interaction. In this energy region, the fluctuations are expected to play a 
decisivly important role for the energy determination of muon neutrino at extremely high energy ( over 
10^^ eV). The fluctuations in this energy region come from (a) the LPM effect due to the muon and (b) 
the LPM effect due to electron (electromagnetic cascade shower). 

In this paper, we would like to discuss the case (a) which has never been discussed before. Without 
exact consideration of the LPM effect in both muon and electron (photon) we could not determine the 
energy of the muon at extremely high energies. 



[l| and in Antarctic 



2 Theory: derivation of cross-section 



At first we remind here the semiclassical picture of LPM. The radiation of a relativistic particle in matter 
develops in a large region of space along its momentum. The characteristic size of this region can be easily 
estimated. When the muon is of a sufficiently high momentum, q the longitudinal momentum carried by 
the virtual photon, becomes very small. 



q=P^.-p'^.-k=^Ef,-ml- ^E'^ - m2 - E^, (1) 

where p^,p'^,E^ and E'^ are the muon momentum and energy before and after the interaction, respec- 
tively, and E^ is the photon energy and 7is Lorentz factor. For high energy muons this simplifies to 



2E^{E^^E^) 2-fl-u 

where u — E^/ E^ is the fractional energy of the radiated photon. This momentum transfer can be very 
small. Then, by the uncertainty principle, the virtual photon exchanges the distance Ic 

h 2El 

lc-- = -TTT (3) 

For example, for a 100 TeV muon emitting a 1 GeV photon, q is 0.56 • 10~^eV/c and the coherence 
length is 3.5 • 10~^ cm. 

The coherence length can also be interpreted as the length in which stripping of a photon from the 
muon which radiates it occurs. 

The LPM effect appears when one considers that the muon must be undisturbed while it traverses 
this distance. One factor that can disturb the muon and disrupt the bremsstrahlung is multiple Coulomb 
scattering. Semiclassically, if the muon multiple scatters by an angle Og, greater than the angle made 
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by the bremsstrahlung photon, 9br ~ l/7i then the bremsstrahlung is suppressed. The average multiple 
scattering angle in a nuclear medium is 

where Eg = ATi/am = 21 MeV is the characteristic energy, I is the scatter thickness, and Xq is the 
radiation length. The LPM effect becomes important when 6s is larger than Ohr- This occurs for 

/X i 

S j V ^0 ^ 7 

If Ic is larger than the size of an atom, it is necessary to take into account the interaction of the 
incident muon not only with the nucleus of the atom but also with the atomic electrons. This interaction 
is taken into account as the screening effect in the Bethe-Heitler theory of bremsstrahlung. If Ic exceeds 
the average distance between the atoms of medium, then it is necessary to take into consideration the 
influence both of the atomic electrons and of many atoms on the particle radiation process. 

The coherence length increases with decrease of the frequency of the radiated photon, and for suffi- 
ciently small it can reach macroscopic dimensions. Therefore, the LPM effect modifies the soft part 

— 1/2 

of bremsstrahlung spectrum. For the low energy photon, the photon spectrum is proportional to E^-^ , 
in contrast to the 1/E^ Bethe-Heitler spectrum. 

An analogous effect occurs for pair creation by a high energy photon. In pair creation, the LPM 
energy threshold is determined by the lepton with the low energy. Because of this, the pair creation 
suppression begins at more higher energies than bremsstrahlung suppression. 

At low bremsstrahlung photon energies dielectric suppression effect occurs (also known as the longi- 
tudinal density effect). This influence can be taken into account by introducing dielectric permeability 
s{E^). Ter-Mikaelian showed, that in the case of soft photons the discount of dielectric permeability also 
leads to decrease of bremsstrahlung probability . Recently these effects have been observed experi- 
mentally 1^ . These phenomena are important in the development of electromagnetic shower accompaning 
the high energy muon. Thus for the modelling of electromagnetic showers it is necessary to use the prob- 
abilities, including LPM effect and the longitudinal density effect. Earlier the discount of LPM effect for 
bremsstrahlung and direct pair production by muon has been carried out in . 

2.1 Bremsstrahlung of muon in a medium 



The appropriate method to consider radiative effects in media is the quasiclassic approach of consid- 
ering the processes of bremsstrahlung and direct pair production Let k = (E^,kj)- be the 
4- momentum of bremsstrahlung photons. 

The probability of radiation of a photon by the high energy muon per a time unit averaged in initial 
muon polarization and summed in final particles polarization has the form 

"^^^J^r^^^j/^Jo '^^«^PHfcWi + T)-x(t))/:(p(< + r),p(t)), (5) 
where C{p{t + T),p{t)) - is a square of matrix elements of muon photon radiation, 

^0 = E^E^/{E^ - E^) - {E^ - P^)/2{E^ - E^), k = E^kJ{E^ - E^) (6) 

This formula determines the radiation probability on the set trajectory and must be averaged in all 
possible trajectories. Averaging is realized due to the distribution function which satisfies to the standard 
kinetic equation 

dF/dt + vdF/df^nv J a{v,v)[F{r,v,t) - F{f,tf ,t)]dv (7) 
with the initial condition 

F{r, r , 0, V, v) = 5{r — r')5{v — v) 
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Here a{v,v') is the scattering cross-section in Couloumb field with the account of full screening 
where x~ is the screening constant in Thomas-Fermi approximation 

In small angle approximation after integration by coordinates the kinetic equation transforms into 
Fokker- Planck equation 

dF/dT + i{a + be^/2)F ^ qAF (9) 

where 



, ^1^1 , E.,E^ 



9i and 02 are minimal and maximal angles of scattering in bremsstrahlung formation process, n is the 
density of scattering centers. The details to solve Fokker-Planck equation can be found in Here we 
show only the final result 

F{e, e',T) = cxp(a(T)^?'2 + l3{T)e9' + 7(t)), (11) 

where 

q;(t) = — ■\/ ib / 8qcoth{y^ 2ibqT) , 
I3{t) = ^/ib/2qsinh-\y/2ib^T), 



7(t) = -iar - ln{sinh{^ 2ibqT)) + d^a{T) + ln{\/ib/8q). 
For C{p ,p) in small angle approximation we have 

c(j^,p)^R^ + R2ee', (12) 

where 



D <^ 1 TJ 



^1) 



El(E^-E^r {E^-E^ 

As a result the quasiclassic approximation calculation gives the probability of bremsstrahlung in the 
following form: 



dW = —q ' 
2n b 



^i?iG(s) + i-i?2$(s) 



3a 



3a 



(13) 



Here s = a/iy/bq, G{s) and $(5) are the famous Migdals functions. 

The indefinity in ln(02/^i) is the result of transition to Fokker-Planck equation in approximation of 
small angles. 

The argument of logarithm contains the indefinite numerical factor. This factor can be defined in the 
limited case (s 1) |^ 

In this case formula (5) goes over to usual bremsstrahlung formula Bethe-Heitler. The comparison of 
the limited expression with the formula of Petrukhin and Shestakov ||l^ , jl^ gives us 



(14) 



where L from llTll is 



L = In 



126m,,Z-2/3/^ 



1 + 18 9Vi(5Z-V3/TO, 

^ Semiclassically, this condition corresponds to 9;,^ ^ 



if Z > 10, 



(15) 
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L = In 



1 + 189ve(5Z-i/3/TO, 



if Z < 10 



(16) 



Here 5 = m?^u/2E^(l — u) is the minimum momentum transfer to the nucleus and e = 2.718. (for the 
electron L = ln(190Z-i/3) ). 

Therefore we can define the value s: 



1/11 



1/2 



27r 7^ (1 — u) nZ'^a^L 



7m 



(17) 



Phenomenological factor ^(s), according to Migdal, equals: 



1, 

l + ln(s)/ln(si) 
2, 



if s > 1 

if si < s < 1 

if s < si 



(18) 



where the value si: 



,1/2 



m, 
m, 



;Zi/Vl89m^, if Z < 10 
;Z2/3/126m^, if Z > 10 



(19) 



The obtained result has the logarithmic accuracy and does not contain Coulomb corrections. The 
use of the more accurate formula |^ as a limit case is difficult because it has not only non-logarithmic 
terms but also corrections proportional to Z and 1/Z. The accuracy of Migdal's calculations is discussed 
at first in |^ . Recently in |2^ the calculation with more accuracy has been performed for electron LPM 
effect. These calculations showed that the Migdal calculations have the 10-15 % accuracy. 

2.2 The longitudinal effect 

Now we consider the longitudinal effect, namely, the influence of media polarization. We can rewrite the 
expression of a from formulae ([ic|): 



mlE^ 



2Ef^ {Efj_ — E. 



2E~f 2E^{E^ — E^) 
1 



^o^tiiEf, - Ey) 



mlE^ 



2E^{Efj, — Ej) 
Therefore we have the additional factor d,, : 



1 



(20) 
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(21) 



This factor changes the expression for bremsstrahlung cross-section (the function <i>(s) <f>(s) /d^, G{s) 
G{s)/d?^ and the value of s ^ • s). 
Finally the cross-section is 



da 
du 



2Z 



i{s)L [lu^G{s) 



d2 



u H — u 

3 3 3 



$(s) 



(22) 



parameter s can be written as 
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, 1 (Elpm u 

We illustrate these both phenomena by some numerical calculations. Parameters of some media are 
collected in Table |l|. 

The condition s < 1 gives us a simple relation for the energy of bremsstrahlung photon, where LPM 
effect is important: 

J^LPM 

The dielectric suppression (effect longitudinal density) is essential at the condition d^> 2: 

< — . 

Thus, for example, the corresponding values in water for muon energy E^ = 10^'^ eV equals -E^ < 25 
MeV (LPM) and E^ < 200 MeV (effect of longitudinal density) 



2.3 Direct electron pair production by muon in a medium 



For the process of direct pair production (DPP) the double differential probability in the total energy 
pair and in the energy of one component (electron or positron) can be obtained, following the paper of 
Ternovskii p^ . The process DPP contains the integration on the virtual photon momentum, therefore 
the answer is expressed through the integral on the Migdal functions. The probability of DPP has the 
following form: 



OTT 



'2\ dE^dE+ 



E'^ 
p 



En 



dW{E^,Ep,E^ 

2 ~1 ( 



El 



El 



En 



where Ep — E^ + E+ and E- are the energies of component pair. 
Let us introduce the following variables: 



P = 



E+ - E^ 
E,+EJ 



E^ 



Ep 



E. 



En 



En 



Thus, we have the following expression for the double differential cross-section: 



Ep 
E,. 



E{s,x) 



(24) 



(25) 



rlz'aHni"-^]^- (26) 



dpdv Sn \ ^1 , 

((1 + /3)(^(s, x) + il + p')B{s, x)) + /3(C7(s, x) + (1 + p^)D{s, x)) + (1 - p^)E{s, x)) 

Here we calculated the Jacobian 

, vdvdp 



and denote 



dE.dE, ^ (27) 



= ^. ^-4°-^^. ™ 

2(1 — v) mg 4(1 — v) 
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The integration limits are: 

^<.<1-^ZV3 (30) 

< IpI < (^1 - v) )^^' 'irnjE.v (31) 

The expressions for coefficients A{s,x), B(s,x),C{s,x), D(s,x)Bu-n.A E{s,x) are given in Appendix A 
0. The argument of logarithm contains the indefinite numerical factor. This factor can be defined in the 
limited case (s ^ 1). In this case formula (7) goes over to usual formula of cross-section direct pair 
production. The comparison of the limited expression with formula of Petruchin and Kokoulin [ pT| , p8| 
gives us 

ln( J) = = L. (32) 



n 



where from fLTl 



Here 



Le « ln(189Z"i/V(l + a:)(l + Ye)) (33) 



Y ^ 5~p'+AP{l + p') 

" 2(1 + 3/3) ln(3 + l/x) - p2 _ 2^(2 - p2) ^ > 

For numerical estimations we can take out the minimal value (1 + x)(l + Yg) 3.25, therefore L = 
ln(614Z^^/'^). The values of Elpm for water, standard rock and lead are listed in Table |[ 

The role of LPM is essential for the process of DPP at 

Elpm 1 



s < 1, i.e.E^ > 



^- 16 v(l-p2)- 

The values of Elpm for bremsstrahlung and for DPP differ from each other only by the quantity L. 
This logarithm is the slowly changeable value depending on the medium properties. 

The energy E"^ is larger for the process of DPP, since the momentum of virtual photons > 4TOg 
for the production of e~^e~ pairs. The LPM is more expressed in the symmetrical case, i.e. p = or 

^ for the direct pair production. 

The probability of production /i^/i^ is 

dW^+^- « -^dWe+e- (35) 
therefore it should not be probably taken into account. 

3 Numerical calculation of the cross-section 
3.1 Bremsstrahlung process 



^ The error from using these expressions gives us 0.7 % accuracy at Efj_ < lO^'^ eV and is less than 15 % accuracy at 
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3.1.1 The longitudinal density effect (Ter-Mikaelian effect) 

The longitudinal density effect ( the LD effect) exists irrespectivly of the absence or presence of the LPM 
effect in the bremsstrahlung process. Therefore, we evaluate the contribution from the LD effect to the 
differential cross section. In Figures 1 to 4, we calculate the differential cross sections for water, standard 
rock, iron and lead, respectively, where BH denotes "without the LPM effect". In these figures, the 
differential cross sections of BH, (BH and LD), LPM and (LPM+BH) are given. It is easily understood 
that the LD effect suppresses much lower energy photon, compared to the LPM effect. And it is also 
understood that the LD effect is more effective in higher density materials. Further it should be noted 
the following: the fact that BH+LD coincide with LPM+LD at smaller u in from Figs.l to 4 meaning 
that the LD effect is only effective in the region at smaller u. However the LD effect is not so effective 
compared to the LPM effect, as in the same situation in the case of electron ^ 

In order to examine the contribution from the LD effect over muon's propagation, we calculate the 
corresponding integral cross section ( the total cross section) due to the LD effect in Figures 5 and 6. 

In Figure 5, we give the ratio of (BH and LD) to BH from 10^ eV to 10^'* eV. In the calculation for 
total cross sections throughout present work, the lower bound for integration is 10^ eV. From the figure 
we could realize that the LD effect is not effective up to 10^^ eV. The LD effect begins to be effective 
around 10^^ eV and makes the integral cross sections decrease to 60 % or more at 10^'' eV in the absence 
of the LPM effect. However, we notice that not a small apparent decrease of the ratio exclusively comes 
from infra-red divergence due to the Bethe-Heitler's cross section which never influences over the real 
behavior of higher energy muon and besides the LD effect is negligible compared to the LPM effect at 
such higher energies , which are shown in Figure 6. Comparing Figure 5 with Figure 6, it is easy to 
understand that (1) the LD effect is only effective and the LPM effect is never effective up to around 
10^^ eV, (2) the competition occurs between the LD effect and the LPM effect from 10^'* eV to 10^" eV, 
and (3) the LD effect is completely negligible at the energies more than 10^" eV. 



3.1.2 The LPM effect 

From Figures I to 4, we understand that the LPM effect is more effective in higher density materials 
and becomes strong at more higher primary energy of muon. The radiated photon at lower energies is 
proportional to ^ 1 / ^JE^ , while that due to BH is proportional to ^ l/ii'7, which is essentially the same 
as in the case of an electron. 

In Figures 7 to 10, we give mean free paths due to the combination of different effects in water, standard 
rock, iron and lead, respectively. The mean free path of an electron due to the LPM effect increases as 
primary enegry of muon increases and as the density of material is higher. Roughly speaking, elongations 
of the mean free paths are 20, 30, 70, 120 times as large as compared to the BH case in water, standard 
rock, iron and lead for the primary energy muon of 10^* eV, respectively. The mean free paths of muons 
with LPM effect are given in Figure 11, while those of the muon with both LPM and LD effects in Figure 
12. Comparison of Figure 11 and Figure 12 shows that the mean free path of muon with both LPM 
effect and the LD effect increases by about fifty percent over 10^^ eV and 10^" eV, which depends on the 
materials, due to the existence of the LD effect compared to the LPM effect only. 



3.2 Direct electron-positron pair production process 

As shown in the Appendix B, the LPM effect in the direct electron pair production process could be 
completely negligible up to 10^"* eV even for lead. Therefore, we calculate the cross sections without 
the LPM effect only. In Figures 13-1 and 13-2, we give the virtual photon energy spectrum for primary 
energies of 10^^ eV and 10^** eV for the materials of lead, iron, standard rock and water, respectively. 
After the integration by the energies of virtual photon, we obtain total cross sections due to the direct 
electron pair productions. In Figure 14, we give the total cross sections for lead, iron, standard rock and 

^See subsequent paper 
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water. In Figure 15, we give the mean free paths of direct pair electrons for the same materials as in 
Figure 14. 



4 Effective energy loss 



Here let us define the effective energy loss of muon in the bremsstrahlung process [2^ by examining the 
influence over "effective energy loss " to primary energy of muon. 

TliEf,) ^ Ef,— J du-u- —a^^^^^^^^{u,Ef,), (36) 
where 4^cr, (u, En) are the differential cross sections with and without LPM effect. Here, we take 

Umin = lO*" eV and Umax = 1 " m^c^ / Ef^ . 

In the effective energy loss by muon, the influence of the longitudinal density effect could be completely 
neglected over 10^ eV to 10^'* eV in both the BH and the LPM cases. 

In Figures 16-1 and 16-2 the effective energy loss due to bremsstrahlung process are given in both 
cases with and without the LPM effect in water and iron, and standard rock and lead, respectively. To 
clear the LPM effect on the effective energy loss, we give the ratio of effective energy loss with the LPM 
effect over that without the LPM effect in the following. 

Umax 

I du-u - -£cr^^^,{u,Ef,) 

riE,) = ^-^^ (37) 

/ du-u - -^cr^„{u,Ef,) 

In Figure 17, the ratio calculated in ( |36| ) is given. It is shown in the figure that the LPM effect becomes 
effective over 10^" eV, depending on the materials. 

Also let us define the average energy transfer to the virtual photon due to direct pair electron pro- 
duction process in the following way. 



^{E^)^E^- dv V -—{v,E^)dp, (38) 



A- Jy^^^ Jp^^^ dpdv ' 

where -^^{v, Ep)dp is a double differential cross section for p and v in the direct pair production. Here 
the lower and upper bounds for integration are defined in (|3^), (|3l|). The calculated results are given in 
Figure 18. As the LPM effect could be neglected completely in the direct electron production process 



shown in the previous section, the numerical results of (37) are given in the BH cases only. 



5 Depth intensity relation of muons at extremely high energies 



Studies on the depth intensity relation of muon at deeper depth in both water and rock are indispensable 
for elucidation of physics related to an extremely high energy muon. Depth intensity relation of muon is 
determined by the processes of bremsstrahlung, direct pair production, nuclear interaction and ionization 
loss. For the studies on the depth intensity relation of muon at extremely high energies, examination of 
the fluctuation effect from the bremsstrahlung process should be carefully made, because it is a strong 
source of the fluctuation. The following calculations are made to examine a fluctuation effect from 
bremsstrahlung. To obtain physically meaningful results at extremely high energies in the calculation 
on the depth intensity relation of muon one should consider not only bremsstrahlung but also the direct 
electron-positron pair production, nuclear interaction and ionization loss down to lower energies, for 
example at 1 GeV. 
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To examine the degree of the fluctuation from bremsstrahhmg process at extremely high energies, we 
calculate the depth intensity relation of muons , taking into account the bremsstrahlung process only and 
keeping Eprim/ E^ut to be 1000, where Eprim denotes the primary energy of the muon and Ecut denotes 
cut off energy of muon in the calculation and the related quantities. 

In Figures 19-1 and 19-2, we give the average value of muon with and without the LPM effect in 
both standard rock and water. We give the average energy of muon with primary energies of lO'^'^eV 
and lG^"^ey. If the LPM effect becomes effective, the cross section decreases so that the average energy 
of muon becomes high. From Figures 19-1 and 19-2, we could understand that [1] the average energy of 
muon in standard rock is higher than that of water, [2] the LPM effect is more effective in lO^^eV than 
in IQ^^eV in both cases with and without the LPM effect. 

In Figure 20-1 and Figure 20-2, survival probabilities of muon are given in rock and water, respectively. 
Conditions given to the calculations are exactly the same as in Figures 19-2 and 21. If the LPM effect 
becomes essential, the cross section decreases so that survival probabilities become large. We could 
understand the following from both Figures . [1] Survival probability in the rock is larger than that in 
water. [2] The LPM effect is much more effective in lO^^eV than in lO^^ey, [3] it gives much higher 
survival probabilities for every depth. 

In both Figures 21-1 and 21-2, range fluctuation of muon are given in both standard rock and water. 
From both Figures, we could conclude the following: comparing the contrast of LPM 24-21 and LPM 
23-20 in Figure 21-1 with those in Figure 21-2, the LPM effect is more essential in standard rock than 
that in water. 

Prom Figure 21-2, Table 3 and Table 4, it is very interesting to note that the average values and their 
standard deviations arc larger in water than in standard rock both with and without the LPM effect, as 
far as the fluctuation of range distributions are concerned. 

In Figures 25-1,25-2 and 25-3, we show energy spectrum in water at the depths lO^g/cm?, 4- lO^g/cm? 
and 8- 10^ g/crn^ for primary energies of lO'^'^eV, lO^^eV^ and lO^^'eF, in the case with the LPM effect and 
without the LPM effect(BH), respectively. From Figure 22-1, we could easily understand that the LPM 
effect is not essential almost over every depth and it shows little effect at depth 8 • 10^ g/cm^, because 
the chances to lose energy increase at larger depths than that at shallower depth so that the LPM effect 
begins to work. From Figure 23-2, we could understand that the LPM effect begins to become effective 
around lO^^eV From 23-3, we understand that the LPM effect becomes effective strongly, particularly at 
the depth 8 • lO^.g/cm^. 

In Figures 23-1,23-2 and 23-3, we show the corresponding spectrum to Figures 22-1,22-2 and 22-3 in 
the standard rock. From the comparison of 23-1, 23-2, and 23-3 with Figures 22-1, 22-2 and 22-3, we 
could understand that the LPM effect is much effective in the standard rock than in water. 



6 General Discussion 



The practical purpose of the study of the muon with LPM effect for neutrino astrophysics is to offer 
theoretical means by which we could determine the energy of astronomical sources with extremely high 
energies . Such extremely high energy sources are expected to be found in the earth through the detection 
of the extremely high energy muon due to muon neutrino. The extremely high energy muon is suffered 
from strong fluctuation due to the LPM effect so that the determination of energy is so difficult. Therefore, 
the detailed studies on extremely high energy muon are requested for the purpose. Another purpose of 
the study on extremely high energy muon may be to clarify the mechanism of charm hadronproduction. 

From the results thus obtained in too much simplified way, we try to conjecture results which cor- 
respond to the real experiments. To perform calculations of the depth intensity relation of muon which 
are really useful for the analysis of the experimental data at extremely high energy, if really exist, we 
must consider the direct electron pair production and nuclear interaction except the bremsstrahlung, and 
energy spectrum of muon at sea level and must make a calculation down to too much lower energy which 
corresponds to the real experiment compared to the present minimum energy lO^^eV, say, lO^eV. Nev- 
ertheless, we could conjecture the following: the LPM effect becomes effective at extremely high energies 
and it appears strong at larger depth due to strong fluctuation effect, because characteristics of relation 
between primary muon with extremely high energy, say lO^^eV and that with still extremely high energy 
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cut off, say lO^^eV, will be surely kept at the relation between primary muon with extremely high en- 
ergy and that with too lower energy , say, lO^eV which corresponds to the real experimental conditions. 
Further, over energy of lO^^eV, the flux of prompt muon at sea level exceeds over that of conventional 
muon which is analyzed in the usual experimental data. Therefore, we could say that the LPM muon is 
appeared in the flux of prompt muon and the depth intensity relation of muon at extremely high energy 
is closely related to the analysis on the mechanism of the charm hadroproduction. 
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Appendix A 



A On calculation of Ternovskii functions 



The integral representations of coeflacients A{s, x), B{s, x), C{s, x), D{s, x), E{s, x) are given by 



, G(sz)dz , $(.sz)dz , r- 

C{s,x)=x — —, D{s,x) = x — —, E{s,x) = 



G{sz)dz 



Now we write the approximation formulae for the functions A(s, x), i?(s, a;), C(s, x), D[s,x),E{s,x). 
The error from using these expressions gives us 0.7 % accuracy at < 10^^ eV and is less than 15 % 
accuracy at > 10^^ eV. Here we used the following simple expression for the functions G{s) and $(s). 



6s + 1' ' ' (6s)2 + 1 

18s2 / 72s2 \ 36s2(l + a;)2 + l Ses^ 
36s2 + l 1^ +36s2 + l7 36^ 36s2 + 1 

216s3 / 36s2 - 1 

+ „„ o - 1 



36s2 + 1 V 36s2 



jxj (^tan(6s(a; + 1)) - I) 



„, , 6s 6sa; \, 6s(l + a;) + 1 6s ^ , „^ 

B(s,x) = IH In— ^— ^— ^ (A -2) 

^ ' ^ 6s + 1 V 6s + 1 y 6sx 6s + 1 ^ ' 

, (36s2)2a; , 36s2(l+:r)2 + l 365^ 
C{s,x) = - .^^. , ^,, ln —5-2 + 



(36s2-|-l)2 36s2a;2 36s2 + 1 

216s3(36s2 - 1) / , ,^ , 7r\ 
- (36s2 + l)2 ' -(Aian{^s{x^X))--^ 

^, , 6s 36s2a; , 6s(l + a;) + 1 
^(^'^) = 6^-(6^^" 6^^ 

E{s, x) = 6s - Atan{%8{x + 1))) 

At s 1/(1 + 0;) in complete agreement with the results of the theory neglecting effects of the medium 
we obtain 
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A{s, x) = {l + 2x) In y^ + lj - 2, B{s, x) = (1 + a;) In + ~ 1' 

C{s,x) = ^-^-^-2x\n(l + -], {A- 3) 

1 + X \ X ) 

i:»(s,j:) = 1 -xln ( 1 + - ) , E{s,x)= ^ 



X I 1 + X 



And at s <C 1/(1 + x) we obtain the expressions 



A(s, x) = — 36s^ ln(6sx), i3(s, = — 6s ln(6sx), C(s, = 36s^, 

D{s,x)=&s, E{s,x) = 3tts. (^-4) 



Appendix B 



Here we show that the LPM effect could be completely negligible for the direct electron pair production 
process up to 10^^ eV even for lead, by using the Ternovskii functions defined in the Appendix A. The 
double differential cross section for the direct electron pair production (DPP) process is given in (p7|). The 
LPM effect is included in the Ternovskii functions from A(s,x) to E{s,x) defined (A-1). The numerical 
orders of all the Ternovskii functions are the same. Therefore we consider the degree of the LPM effect 
in the function A{s, x) in the Appendix A, as an example. 

A{s,x) is given as follows: 

Here, we write down 

G{s) = l-{l-G[s)} [B-2) 
If we insert (B-2) into (B-1), then we obtain 

/"^ {z-x~ l)dz f°° {z - X - - G{sz)]dz 

The first term of the right side (B-3) represents the term without LPM effect, the second term in (B-3) 
describes the influence of the LPM effect. At G{sz) = 1 the second term vanishes. Therefore, we could 
rewrite the (B-1) in the following way: 

A{s,x) ^ A{BH) + A{LPM), {B - A) 

where 

Jl+x Z\Z-IY 

AiLPM) ^ - r i^-^-m-Gisz)]dz 

Jl+x z2(z-l)2 

As (1 — G{sz)) in (B-6) is a non-negative monotonically decreasing function, A{LPM) is a non-positive 
function, which decreases the cross section as the logical consequence of the LPM effect. If the absolute 
value of the A{LPM) is negligible compared to A{BH) , then, it is equivalent to the neglect of the LPM 
effect. 

Now, let us examine the degree of the LPM effect , calculating the A{s, x),B{s, x), C{s, x), D{s, x),E{s, x) 
for Ef^ ~ 10^^ eV in the lead. Thus, as the LPM effect is the strongest in the lead among four kinds of 
materials under consideration, it is enough to consider the case of the lead. 

In Table B-1 we give A{BH),A{LPM) - E{BH),E{LPM). As it is clearly seen from the Table B-1, 
all the functions A{LPM) — E{LPM) are completely negligible for A{BH) — E{BH), respectively over 
all V and p. This is the reason why we assert that the LPM effect could be completely negligible for the 
direct pair production process up to 10^"* eV. 

In Table B-1, we give numerical values of A{BH), A{LPM) for different v and p in the case of 
lead with As you understand from these numerical values, the LPM effect could be completely 

negligible even in the lead with 
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Figures and Table Captions 



Figure 1 Differential cross sections for brcmsstrahlung processes in water. LPM denotes the case with 
the LPM effect and LPM+LD denotes the case with longitudinal effect (LD) in addition to the LPM 
effect, while BH denotes the case without the LPM effect (the Bethe-Heitler case) and BH+LPM 
effect denotes the case without the LPM effect and with the LD effect, [a] denotes lO-'^^eV, the 
primary energy of the muon and [b] to [ji] denote 10^^ to lO^^eV, primary energy of the muons, 
respectively. 

Figure 2 Differential cross sections for bremsstrahlung processes in standard rock. The notations in the 
figure are the same as in Figure 1 

Figure 3 Differential cross sections for bremsstrahlung processes in iron. The notations in the figure are 
the same as in Figure 1. 

Figure 4 Differential cross sections for bremsstrahlung processes in lead The notations in the figure are 

the same as in Figure 1. 

Figure 5 Ratio of the total cross sections (Tbh+ld to ubh for various materials. 
Figure 6 Ratio of the total cross sections (Jlpm+ld to (Tlpm for various materials. 

Figure 7 Mean free paths for bremsstrahlung processes in water The notations in the figure are the 

same as in Figure 1. 

Figure 8 Mean free paths for bremsstrahlung processes in standard rock The notations in the figure are 
the same as in Figure 1. 

Figure 9 Mean free paths for bremsstrahlung processes in iron. The notations in the figure are the same 
as in Figure 1. 

Figure 10 Mean free paths for bremsstrahlung processes in lead. The notations in the figure are the 
same as in Figure 1. 

Figure 11 Mean free paths for bremsstrahlung processes with the LPM effect for various materials. 

Figure 12 Mean free paths for the bremsstrahlung process with the LPM and the LD effect for various 
materials. 

Figure 13-1 Differential virtual gamma rays energy spectrum for direct electron pair production pro- 
cesses for the primary energy with lO^^eV 

Figure 13-2 Differential virtual gamma rays energy spectrum for direct electron pair production pro- 
cesses for the primary energy with lO'^'^eV. 

Figure 14 Total cross sections for direct electron pair production process for various materials 

Figure 15 Mean free paths for direct electron pair production process for various materials. 

Figure 16-1 Effective energy loss for bremsstrahlung process due to both with and without LPM effct 
for the iron and water. 

Figure 16-2 Effective energy loss for bremsstrahlung process due to both with and without LPM effct 
for the lead and s.r. 

Figure 17 The Ratio of the effective energy loss with the LPM effect to without the LPM effect for 
bremsstrahlung process for various materials 

Figure 18 Effective energy loss for direct electron pair production process for various materials. 
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Figure 19-1 Decrease of average energy of muon in the standard rock Ordinate axis means average 
energy of muon divided the sampUng number of muons for both the cases with and without the 
LPM effect. LPM 24-21 means primary energy of muon with lO'^^eV and the cut off energy of 
muon with lO'^^eV with the LPM effect, while BH 24-21 correspond to the case without the LPM 
effect (Bethe-Heitler case). LPM 23-29 means primary energy of muon with lO^^eV and the cut off 
energy of muon with lO^^eV with the LPM effect , while BH 23-20- means the primary energy of 
muon and the cut off energy of muon with lO'^^eV without the LPM effect . 

Figure 19-2 Decrease of average energy of muon in water. The date correspond to those in Figure 19-1. 

Figure 20-1 Survival probability of muon in standard rock LPM 24-21 and other notations are the same 
as in Figure 19-1. 

Figure 20-2 Survival probability of muon in water. The date correspond to those in Figure 19-1. 

Figure 21-1 Range distribution of muon in standard rock LPM 24-21 and other notations are the same 

as in Figure 19-1 

Figure 21-2 Range distribution of muon in water. Notations in the figure are the same as in Figure 
19-1. The area covered with the range distribution is normalized to unity. 

Figure 22-1 Differential energy spectrum in water : LPM 22-19 and BH 22-19 are of the same meaning 
as in Figure 19-1. 

Figure 22-2 Differential energy spectrum in water : LPM 23-20 and BH 23-20 are of the same meaning 
as in Figure 19-1. 

Figure 22-3 Differential energy spectrum in water : LPM 24-21 and BH 24-21 are of the same meaning 
as in Figure 19-1. 

Figure 23-1 Differential energy spectrum in standard rock : LPM 22-19 and BH 22-19 are of the same 

meaning as in Figure 19-1. 

Figure 23-2 Differential energy spectrum in standard rock : LPM 23-20 and BH 23-20 are of the same 
meaning as in Figure 19-1. 

Figure 23-3 Differential energy spectrum in standard rock : LPM 24-21 and BH 24-21 are of the same 
meaning as in Figure 19-1. 

Table 1 Some parameters for the bremsstrahlung process 

Table 2 Some parameters for the direct electron pair production process 

Table 3 Mean values and their standard deviations of range fluctuation of muons in water. 

Table 4 Moan values and their standard deviations of range fluctuation of muons in standard rock. 
Table B-1 A{BH) and A{LPM) in the Ternovskn functions. 
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Table 1: 





Water {H^O) 


St. rock 


Lead (P6) 


Fe 


z 


7.23 


11 


82 


26 


A 


14.3 


22 


207.2 


55.85 


p{g ■ cm^'^) 


1.0 


2.65 


11.34 


7.86 




4.211 • lO""^ 


7.254- 10^^ 


3.296 • 10^^ 


8.475 • 10^^ 


L 


8.85 


8.57 


7.23 


7.99 


ELPM{eV) 


2.49 • lO^'^ 


6.44 • 10^^ 


3.02 • 10^^ 


1.06 • 10^^ 




20.49 


33.18 


61.06 


55.14 



Table 2: 





Water {H2O) 


St. rock 


Lead (Pb) 


Iron (Fe) 


L 


5.76 


5.62 


4.95 


5.33 


ELPMieV) 


3.84 • 10^4 


0.98 • 10^* 


4.41 • 10^^ 


1.59 • 10^^ 



Table 3: 



Muon energy 


BH 


LPM 


Eprim[eV] Ecut[eV] 


< i? > a 


< R> a 


10^* 10'^' 


34.28 14.60 


38.44 16.93 




32.71 14.53 


33.79 15.00 


1019 


33.85 14.97 


33.95 15.04 


10^1 ioi« 


35.00 15.39 


34.99 15.42 


10^0 ^Ql7 


36.09 15.77 


36.07 15.80 



Tal)l(^ 4: 



Muon energy 


BH 


LPM 


Eprim[eV] Ecut[eV] 


< R> a 


< i? > cr 


10^4 10^^ 


26.36 11.23 


35.30 15.82 




25.15 11.17 


27.32 12.09 


10^^ lOiy 


26.03 11.51 


26.31 11.66 


10^1 ioi« 


26.91 11.84 


26.91 11.86 


10^0 ^Ql7 


27.75 12.13 


27.74 12.15 
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Table 5: 



V 


p 


A (BH) 


A (LPM) 


1.022 


10" 


18 


0.000 


lO'J 




7 


157263 


10 


0.000000 • 10" 


1.000 


10" 


17 


0.000 


10° 




6 


701101 


10 


0.000000 • 10" 


1.000 


10" 


17 


5.000 


10" 


1 


6 


729869 


10 


0.000000 • 10° 


1.000 


10" 


17 


9.475 


10- 


1 


6 


929182 


10 


0.000000- 10° 


1.000 


10" 


16 


0.000 


10° 




6 


240584 


10 


0.000000 • 10° 


1.000 


10" 


16 


5.000 


10" 


1 


6 


269352 


10 


0.000000 • 10° 


1.000 


10" 


16 


9.949 


10- 


1 


6 


698923 


10 


0.000000- 10° 


1.000 


10" 


15 


0.000 


10° 




5 


780067 


10 


0.000000 • 10° 


1.000 


10" 


15 


5.000 


10" 


1 


5 


808835 


10 


0.000000 • 10° 


1.000 


10~ 


15 


9.995 


10- 


1 


6 


468665 


10 


0.000000- 10° 


1.000 


10" 


14 


0.000 


10° 




5 


319550 


10 


0.000000 • 10° 


1.000 


10" 


14 


5.000 


10" 


1 


5 


348318 


10 


0.000000 • 10° 


1.000 


10- 


14 


9.999 


10- 


1 


6 


238406 


10 


0.000000- 10° 


1.000 


10" 


13 


0.000 


10° 




4 


859033 


10 


0.000000 • 10° 


1.000 


10" 


13 


5.000 


10- 


1 


4 


887801 


10 


0.000000 • 10° 


1.000 


10- 


13 


1.000 


10" 




6 


008148 


10 


0.000000- 10" 


1.000 


10" 


12 


0.000 


10° 




4 


398515 


10 


0.000000 • 10" 


1.000 


10" 


1 o 


5.000 


10" 


1 


4 


427284 


10 


0.000000 • 10° 


1.000 


10- 


12 


1.000 


10" 




5 


777889 


10 


0.000000 - 10" 


1.000 


10" 


11 


0.000 


10° 




3 


937998 


10 


0.000000 • 10" 


1.000 


10" 


11 


5.000 


10" 


1 


3 


966767 


10 


0.000000 • 10° 


1.000 


10- 


1 1 


1.000 


10" 




5 


547631 


10 


0.000000- 10° 


1.000 


10" 


10 


0.000 


10" 




3 


477481 


10 


0.000000 • 10° 


1.000 


10" 


10 


5.000 


10" 


1 


3 


506250 


10 


0.000000 • 10° 


1.000 


10- 


10 


1.000 


10" 




5 


317372 


10 


0.000000- 10° 


1.000 


10" 


9 


0.000 


10" 




3 


016964 


10 


0.000000 • 10° 


1.000 


10" 


9 


5.000 


10" 


1 


3 


045733 


10 


0.000000 • 10° 


1.000 


10- 


9 


1.000 


10" 




5 


087114 


10 


0.000000- 10" 


1.000 


10" 


8 


0.000 


10" 




2 


556447 


10 


0.000000 • 10" 


1.000 


10" 


8 


5.000 


10" 


1 


2 


585216 


10 


0.000000 • 10° 


1.000 


10" 


8 


1.000 


10" 




4 


856855 


10 


0.000000- 10" 


1.000 


10" 


7 


0.000 


10" 




2 


095930 


10 


0.000000 • 10" 


1.000 


10" 


7 


5.000 


10" 


1 


2 


124699 


10 


0.000000 • 10° 


1.000 


10" 


7 


1.000 


10" 




4 


626597 


10 


0.000000- 10° 


1.000 


10" 


6 


0.000 


10" 




1 


635413 


10 


0.000000 • 10° 


1.000 


10" 


6 


5.000 


10" 


1 


1 


664182 


10 


0.000000 • 10° 


1.000 


10" 


6 


1.000 


10" 




4 


396338 


10 


0.000000- 10° 


1.000 


10" 


5 


0.000 


10" 




1 


174898 


10 


-3.030322- lO-*" 


1.000 


10" 





5.000 


10" 


± 


1 


203666 


10 


-1.751002- 10"^ 


1.000 


10" 


5 


1.000 


10" 




4 


166079 


10 


0.000000- 10° 


1.000 


10" 


4 


0.000 


10" 




7 


145755 


10" 


-1.698723- 10-'' 


1.000 


10" 


4 


5.000 


10" 


1 


7 


432968 


10° 


-1.002189- 10"^ 


1.000 


10" 


4 


1.000 


10" 




3 


935811 


10 


0.000000- 10° 


1.000 


10^ 




0.000 


10" 




2 


645587 


10" 


-4.427908- 10^ 


1.000 


10^ 




5.000 


10" 


1 


2 


910863 


10° 


-2.918965- 10^ 


1.000 


10^ 




1.000 


10" 




3 


705456 


10 


0.000000- 10° 


1.000 


10^ 




0.000 


10" 




7 


112332 


10^ 


-2.047297- 10-"* 


1.000 


10 




5.000 


10" 


1 


1 


066582 


1 n — 1 

10 


-3.217814- 10"^ 


1.000 


102 




1.000 


10" 




3 


474261 


10 


U.UUUUUU - iU 


1.000 


10- 


1 


0.000 


10" 




1 


171890 


10-5 


-2.204098- 10-13 


1.000 


10- 


1 


5.000 


10- 


1 


2 


077567 


10-5 


-7.148156- 10-13 


1.000 


10- 


1 


1.000 


10" 




3 


233732 


10 


0.000000-10° 


1.000 


10° 




0.000 


10" 







000000 


10° 


0.000000- 10° 


1.000 


10° 




5.000 


10- 


1 





000000 


10° 


0.000000 • 10° 


1.000 


10° 




1.000 


10° 




2 


363193 


10-1° 


0.000000 • 10° 
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BH24-21 

BH23-20 

LPM24-21 

LPM23-20 

S.R. 



50 100 [x10n 

Depth [g*cnn~^] 




Ea [eV] 




Depth [g*cnn ^] 



0.06- 



0.04- 



0.02- 



0- 



L 







b- 



_ I 



J I \ I I L 



50 



-2n 



Depth [g*cnn ] 



1 1 — 

BH24-21 
BH23-20 " 

LPM24-21 

- LPM23-20 

Water 




100 



[x10^^] 




10 



12 



10 



16 



10 



20 



10 



24 



